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A mathematical model is developed that captures the transport of liquid wa- 
ter in hardened concrete, as well as the chemical reactions that occur between 
the imbibed water and the residual calcium silicate compounds residing in the 
^ ' porous concrete matrix. The main hypothesis in this model is that the reac- 

tion product - calcium silicate hydrate gel - clogs the pores within the concrete 
thereby hindering water transport. Numerical simulations are employed to de- 
I termine the sensitivity of the model solution to changes in various physical 

' parameters, and compare to experimental results available in the literature. 

cn ■ 1 Introduction 

. Concrete is a ubiquitous construction material that derives its utility from a 

QQ ' combination of strength, versatility and relatively low cost. In fact, concrete 

. is the second most consumed material on the planet next to water [?7]. The 

primary ingredients that go into the making of concrete are water, Portland 
cement (a fine powder consisting primarily of calcium silicate compounds) , and 
^ I solid aggregates such as sand and gravel. When mixed together, these ingredi- 

^ . ents undergo a complex physico-chemical transformation which can be divided 

into a number of discrete steps: an initial hydration stage that occurs over a 
period of hours or days; a drying/curing period that can require months or even 
years to complete; and additional reactions arising from carbonation and various 
degradation processes that typically also occur over very long time periods. 

Mathematical modelling of transport and reaction in concrete has been the 
subject of a large number of papers in the scientific and engineering literature. 
The earliest study that we are aware of which treats water transport in concrete 
as a nonlinear diffusion process is that of Bazant and Najjar 6]. Later work 
considered the additional effect of transport and reaction of chemical species in 
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the context of initial cement hydration [T31 SIl [S] or concrete carbonation [351 
[37l [T9] . We remark that most of these models assume a constant porosity even 
though experimental evidence overwhelmingly suggests that the pore structure 
varies significantly over time owing to reactant consumption, crystallization, and 
swelling of products throughout the various stages of concrete hydration [531 
H3] . In fact, the only models we are aware of that allow for a variable porosity 
are in the context of carbonation, where Meier et al. [37] specify the porosity as 
a given decaying exponential function of time, while Bary and Sellier |5] allow 
the porosity to depend on the solution via changes in the pore volume from 
solidified reaction products. 

We focus here on a later stage in the life of concrete, namely the process 
of re-wetting or "secondary hydration" in which hardened and cured concrete 
experiences imbibition of water, due to periodic rainfall for example. The pro- 
portion of reactive silicates in the cement that are consumed during the initial 
hydration reactions (called the degree of hydration) is typically on the order of 
50% |5] ; consequently, there are significant levels of residual reactants remaining 
in hardened concrete and so the effect of secondary reactions occurring during 
re-wetting cannot be ignored. This study is motivated by the experimental work 
of Barrita et al. [H [3] and Taylor et al. [33] who placed dry concrete samples in a 
liquid bath and carefully observed the progress of the subsequent wetting front. 
They found that when a non-reactive liquid such as isopropanol was used, the 
front speed was proportional to the square root of time as predicted by nonlin- 
ear diffusion analysis. When water was used instead, the wetting front moved 
more slowly than the theory predicted and in some cases stalled completely - an 
effect that is usually referred to as anomalous diffusion. Hall [24j has suggested 
that this effect is due to physico-chemical interactions between the wetting fluid 
and the porous solid. It is natural to hypothesize therefore that the reduction in 
wetting front speed arises from residual calcium silicates in the porous matrix 
reacting with water to form calcium silicate hydrate or C-S-H, which precipi- 
tates in the form of a gel that clogs the pores in the concrete; this hypothesis is 
supported by the results in 

Observations of anomalous diffusion have been reported in [35] where the 
authors proposed instead that deviations in wetting front speed can be mod- 
eled using a modified (non-Darcian) porous transport equation. This approach 
provides a reasonable match with experiments and gives rise to a new and poten- 
tially interesting class of nonlinear diffusion equations and scaling laws; however, 
there is no direct support for this model in terms of a physical mechanism for 
concrete hydration. In a related study |35] , another model is proposed that 
includes an explicit time-dependence in the water diffusion coefficient. They 
showed that by assuming the cumulative deposition follows a power law in time, 
they could reproduce similar clogging results; unfortunately, it is not at all clear 
how one would obtain the power law coefficients in a given wetting scenario. 

Some authors have addressed clogging phenomena in the context of concrete 
carbonation, such as Saetta et al. |42| who incorporated a functional depen- 
dence on the carbonate concentration into the transport coefficients of their 
model. Meier et al. [37j also employed an empirical approach, but instead they 
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assumed the porosity decays as a given exponential function of time, which has 
the disadvantage that there is no direct coupHng between water transport and 
the precipitated reaction products that are causing the actual clogging. Re- 
lated work on self-desiccation (or internal drying) during initial hydration and 
its connection with autogenous shrinkage have been studied using pore-level 
microstructure simulations [10. 

In this paper, we develop a model that aims to test the hypothesis that in- 
corporating the chemistry of residual cement constituents and the effect of the 
resulting C-S-H gel formation on pore structure can explain the apparent clog- 
ging effects observed in concrete re- wetting experiments. We begin in Section [5] 
with a brief overview of cement chemistry and the physico-chemical changes 
that occur in cement during hydration. We develop the mathematical model in 
Section [3] using a macroscopic approach that is motivated by the clogging mod- 
els developed in the bioremediation literature (see for example |16j ) wherein 
the accumulation of biomass - analogous to cement hydration products - is 
responsible for the reduction in porosity. Numerical simulations of the result- 
ing system of nonlinear partial differential equations are performed in Section!?] 
and the results are compared with experiments. We show that our model cap- 
tures observed clogging behaviour both qualitatively and quantitatively with a 
minimum of parameter fitting, and we explain in Section [5] how these results 
might be generalized in future to handle a range of other related phenomena in 
concrete transport. 

2 Overview of cement chemistry 

While this paper is not concerned directly with the primary hydration of cement, 
the same hydration reactions occur during the re-wetting phase when residual 
unhydrated silicates remaining in the hardened concrete matrix are exposed to 
water. We will therefore begin by presenting some background information on 
the process of cement hydration that is drawn largely from [33l [11] . Portland 
cement is the key binding agent in concrete and has as its major constituents 
tricalcium silicate (3CaO • Si02, commonly referred to as alite) and dicalcium 
silicate (2CaO ■ Si02, known as belite) which make up approximately 50% and 
25% respectively of dry cement by mass. The remaining 25% consists primarily 
of tricalcium aluminate, tetracalcium aluminoferrite, and gypsum, with smaller 
amounts of certain other admixtures whose purpose is to influence such proper- 
ties as strength, flexibility, setting time, etc. In this paper we will concentrate 
solely on the two primary constituents, alite and belite. 

The cement powder is mixed with water and aggregates (sand, gravel and 
crushed stone) to make a workable paste that can then be easily poured or 
molded and left to harden. During the initial hydration stage, the silicates dis- 
solve in and react with the water to form calcium hydroxide or Ca(0H)2, and 
calcium silicate hydrate or C-S-H; the latter notation does not denote a specific 
chemical compound but rather represents a whole family of hydrates having 
Ca/Si ratios that range between 0.6 and 2.0. A significant amount of heat is 
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released during the conversion of alite and belite into C-S-H since the hydra- 
tion reactions are exothermic. Calcium hydroxide and C-S-H precipitate out of 
solution in crystalline form, and these solid precipitates then act as nucleation 
sites that further enhance formation of C-S-H. It is the crystalline or gel form 
of C-S-H that is ultimately responsible for the strength of concrete. 

The hydration process can proceed via several possible reactions, but we 
restrict ourselves to a particular reaction sequence that is employed in both [41j 
and [46j . The mechanism for alite hydration begins with a dissolution phase 

3CaO ■ Si02 + 3H2O ^ 3Ca2+ + 40H- + H2Si04", 
followed by a reaction in solution to form aqueous C-S-H 

HaSiO^" + ^Ca2+ + OH" ^ C-S-H (aq), 
and precipitation of calcium hydroxide according to 

Ca?+ + 20H- ^ Ca(0H)2 

In each chemical formula we have indicated the rate of the reaction by [day^^] 
for i = 1,2,3. 

For the remainder of this paper, we adopt the cement chemistry convention in 
which the following abbreviations are used: C = CaO, S = Si02 and H = H2O. 
Then the overall reaction, leaving out intermediate ionic species, can be written 
in terms of the single formula 

2C3S + 6H ^ C-S-H (aq) + 3CH, (1) 

where Tq represents an overall rate constant for alite. Motivated by the models 
developed in [38l [371 ED i we consider only the simplified kinetics represented 
by (P). We also take the chemical form of C-S-H to be that of C3S2H3, that 
because of the amorphous nature of C-S-H can only be true in some averaged 
sense. A similar formula holds for belite 

2C2S + 4H ^ C-S-H (aq) + CH, (2) 

where we note that ^ [131 IMl E] . Alite is also mainly responsible for 
the early stage strength of the concrete (through approximately the first seven 
days) while belite contributes to the later strength. 

Following [46], we include a precipitation (or deposition) step in which the 
aqueous C-S-H product precipitates out of solution to bind with the porous 
matrix: 

C-S-H (aq) C-S-H (gel), (3) 

where the rate of precipitation is denoted by fcprec [day^^]. We allow for a disso- 
lution process with rate constant k^iss , although in most of our later simulations 



4 



we restrict ourselves to kdiss = so as to be consistent with other models that 
disregard the effect of C-S-H dissolution. 

The hydration chemistry of other cement constituents such as aluminates, 
ferrites, etc. are not considered here because they do not contribute appreciably 
to the porous structure of the concrete [131 [H] ■ Instead, we focus on the effect 
of C-S-H gel on decreasing porosity and hindering moisture transport within 
the porous concrete matrix. 

3 Mathematical model 

We begin by providing a list of primary simplifying assumptions that will permit 
us to reduce the complexity of the governing equations while at the same time 
retaining the essential aspects of the underlying physical and chemical processes: 

1. The length scales under consideration are large enough that the solid 
concrete matrix can be treated as a continuum. Consequently, volume 
fractions and constituent concentrations can be expressed as continuous 
functions of space and time, and the liquid transport is assumed to obey 
Darcy's law. 

2. The concrete sample is long and thin so that transport can be assumed to 
be one-dimensional. This is consistent with many experiments involving 
concrete or other building materials [3] in which the sample under study 
takes the shape of a long cylinder as pictured in Fig. [T^. 

3. Temperature variations and heat of reaction can be ignored. This is a 
reasonable assumption in re-wetting of hardened concrete for which the 
quantities of silicate reactants are much smaller than during the initial 
hydration stage [26] . 

4. Water transport is dominated by capillary action and so gravitational 
effects can be ignored. This assumption is justified by the very small pore 
dimensions that lead to a small Bond number for concrete [35j . 

5. We neglect the dynamics of individual ionic species, which is consistent 
with the work of Bentz llj and others. Nonetheless, we do consider 
separate aqueous and solid phases of C-S-H and include a simple dynamic 
mechanism for precipitation and dissolution which is shown in [46' to be an 
important effect. This choice is motivated by the recognized complexity of 
the C-S-H precipitation/crystallization process [52], that is largely ignored 
by other models of hydration. 

6. The effects of chemical shrinkage and subsequent self-desiccation can be 
significant during initial cement hydration, particularly for high perfor- 
mance concrete [10]. However, we neglect both effects since the samples 
under consideration are normal strength concrete, and residual silicate 
concentrations are much smaller during re-wetting. 



5 



7. Reaction kinetics take a simple form in which the rate has a power-law 
dependence on reactant concentration - a common assumption employed 
in other models [351 ISZ] . 




a. Cylindrical concrete sample, de- b. A zoomed-in view at the pore 
picting the location of a typical wet- scale, showing a representative ele- 
ting front. mentary volume f2. 

Figure 1: Geometry of the ID moisture transport problem. 



3.1 Definition of volume fractions and concentrations 

Consider an elementary volume fl [cm^] pictured in Fig. [T|3 which is divided into 
sub- volumes occupied by the various components of the porous matrix, namely 
the non-gel solids with volume fi^, the precipitated C-S-H gel Qg, liquid water 
flwi and gas/vapour component fly. The pore volume available for transport is 
denoted by flp — Qw + and so the total volume can be written as 

We next define the various volume fractions beginning with the pore volume 
fraction e = flp/fl [cm^ /cm^], which is also known as porosity. The initial poros- 
ity in the absence of C-S-H is denoted by the constant e° = £\t=o = (fl — ^ls)/^. 
The gel volume fraction is Eg = Qg/Q and volumetric water content is 9 = O^j/fi. 
In practice, 9 must satisfy < 9min ^ 9 < 9,nax £°i where 9min is the im- 
mobile or residual water content and 9,nax is the maximum or fully saturated 
value (representing the point beyond which water can no longer penetrate the 
smallest pores). In the context of cement hydration, there are three forms of 
water present: chemically bound, physically bound, and capillary water. The 
quantity 9min corresponds to both physically bound (absorbed) and capillary 
water, which together represent the "evaporable water" that can be removed 
only by forced drying. All volume fractions will in general be functions of both 
position and time owing to variations in the gel, liquid, and gas concentrations. 
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We next define the concentrations of tlie various constituents (in units of 
[g/cm'^]), which are measured relative to the total mass of concrete following 
Ref. [23: 

Ca{x^t) - concentration of C3S in concrete, 
Ci3{x,t) - concentration of C2S in concrete, 
Cq{x,t) - concentration of C-S-H in liquid, 
Cg{x, t) - concentration of sohd C-S-H gel — pgV,g/il — PgSg- 

All solution components are taken to be functions of time t and axial distance 
X, where i > and < x < L. The gel-modified porosity e{x, t) is related to 
C-S-H gel concentration via 

_ ~ Q,S — ilg ^ ^ - _ ^ ^ O _ ^ fA\ 

n n n pg' ^ ' 

where pg is the density of C-S-H in gel form [g/ 
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cm 

3.2 Derivation of the governing equations 

We next derive the differential equations governing the water content 9 and each 
of the constituent concentrations Cq, C/3, Cq and Cg. Conservation of liquid in 
the pores requires that 

f = -|^-.(^-^.)+^^, (5) 
at dx pw rricsh 

where u is the water velocity [cm/ day] and pj^ is its density [g/cm^]. The 
reaction term is scaled by a factor {9 — 9r)^ — max(0 — 6*^,0) which reflects 
the assumption that reactions proceed only when water content is above some 
minimum value 9r- A similar approach was used in studies of concrete carbon- 
ation [371 112] wherein the value of 6*^ is obtained experimentally; in the absence 
of experimental data for hydration, we simply take 9r = 9min ■ The expression, 
Tcsh representing the rate of generation of C-S-H [g/cm'^ day] must be scaled 
here by the ratio of the molar masses of water and C-S-H, rriu] /rricsh- We also 
multiply the reaction term by a stoichiometric coefficient v, which is taken equal 
to 5 so as to balance the rate of generation of water with the averaged rate co- 
efficient for C-S-H coming from Eqs. ((1]) and This and other reaction terms 
are specified later in Section 13.41 

We assume the liquid velocity can be expressed as 

ri9 

^=-D{e,e)-, (6) 

which follows from a simple application of Darcy 's lawQ, where the effective 
diffusivity D{9,e) [cvn? / day] is a function of both water content and porosity 

* Darcy 's law states that the velocity u = K dh/dx, whore h is the pressure head and K is 
the hydraulic conductivity of the medium. Conductivity is known to depend on porosity, and 
for variably saturated media both K and h are typically taken to be functions of saturation, 
such as in the van Genuchten or Brooks-Corey models [2- Therefore, Darcy's law takes the 
form of Eq. ^ with D = K{9, e) dh/de. 
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(for which a specific functional form will be presented in Section [3.5^ . After 
substituting ([6]) into (O, we obtain 



d0 
di 



d_ 

dx 



'^w '^csh 
Pw '^csh 



(7) 



Dissolved alite and belite are advected with the pore liquid as well as be- 
ing affected by diffusion and reaction, and so the corresponding conservation 
equations are 



dt 
dt 



d_ 

dx 
d_ 

dx 



dx 
dCp 
dx 



dx 
dx 



-{9- Br)+ r„ 



(8) 
(9) 



where Dj^ j = a, (3 is the diffusivity of each dissolved constituent. Transport of 
aqueous C-S-H is governed by 



d jdCg) 

dt 



d_ 

dx 



d{uCg 



I") ij' csh 



dx 

kprecCq 



(10) 



5 C, 



9) ' 



where kprec and k^iss are the rates of C-S-H precipitation and dissolution respec- 
tively. The solid C-S-H phase is not affected by advective or diffusive transport 
and so obeys a simple ODE 



dt 



r)^ (fcj 



tree Cq 



(11) 



In summary, the governing equations consist of ([7]), ([8))- pT|) . which enforce 
conservation of water, aqueous species, and solid C-S-H gel, supplemented with 
the relationships ^ and 



3.2.1 Analogy with bioremediation models 

Before presenting the remaining details, it is worthwhile mentioning that there 
is a great deal of similarity between our model for reactive transport in concrete 
and those developed for biofilm growth and bioremediation in the soil sciences 
literature (for example, [121 HSl HHj ■ In the case of bioremediation, bacteria are 
employed in porous aquifers in order to break down some targeted contami- 
nant. Nutrients (typically nitrates) are injected into the ground to activate the 
decontamination process and soil scientists are interested in understanding how 
to encourage the growth of the bacteria in a controlled manner so as to avoid 
clogging the pores in the rock or soil matrix while at the same time maximiz- 
ing the breakdown of contaminant. The governing equations for both problems 
therefore have a similar structure, with a few key differences that we summarize 
below: 
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• In biofilms water is an inert phase, whereas in concrete it participates in 
the reaction. 

• Biological organisms are typically modelled using Monod reaction terms, 
whereas we use power-law kinetics. 

• Biofilms are composed of living cells and so give rise to additional terms 
that encompass cell division and death processes. 

• The microstructurc of biofilms and C-S-H are quite different, but our use 
of a continuum approach means that we can ignore such details. We 
do nonetheless employ the same power-law form (j20p of the permeability 
correction as that used in biofilms. 



3.3 Initial and boundary conditions 

We assume that the concrete sample at the beginning of an experiment is homo- 
geneous in composition and uniformly hydrated so that the initial water content 
and concentrations for < x < L are 

C,(x,0) = C°, Cg{x,0) = C;, 

where the zero superscript denotes a constant initial value. The first condition 
on water content states that the concrete is initially at the minimum value, 
corresponding for example to a sample that is equilibrated in a humidified en- 
vironment but not force-dried. It is reasonable to take the initial C-S-H con- 
centrations C° = Cg = 0, but the alite and belite concentrations are key model 
parameters that depend on the composition of the initial cement mixture. In 
particular, Papadakis et al. [55] calculate the initial concentrations as 

q = (13a) 



where 



Crnix — "TZ J } (-^^^) 

/ c Pcemj Pw ^" Ra/c Pceraj Pagg ~t- 1 



-'cem 



represents the initial concentration of cement before onset of hydration, pc 
is the original cement density, pagg is the particle density of aggregates, Rw/c 
and Ra/c ^-re initial water-to-cement and aggregate-to-cement ratios by mass, 
and ujj is the mass fraction for each constituent j — a, (3. We have modified 
Papadakis et al.'s formula slightly to include the extra factors (1 — fj) where 
fj e [0, 1] represents the fractional degree of hydration of each constituent at 
the end of the hydration/curing stages. 

The cement mixtures investigated in [3] contain significant levels of trical- 
cium aluminate (or C3A, short for SCaO • AI2O3) and no tetracalcium alumi- 
noferrite. Consequently, for the purposes of calculating initial porosity, we also 
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include the effect of C3A, whose initial hydration products further reduce the 
pore space available for transport. Letting and uj^ refer to the mass and 
hydration fractions for C3A, we are led to the following expression for initial 
porosity i38:, Tab. 2]: 

£° = CmixRtu/c/Pw - Cmix (fa^^a AV^ + fpUJp AVp + f^UJ^ /W^) , (14) 

where the first term represents the porosity before onset of hydration and the 
remaining terms encompass the reduction in pore volume owing to hydration 
through parameters AVq = 0.233, /Wp = 0.228 and AF^ = 0.555 (units of 



cm 



l9])- 



We note in passing that the strength of the resulting hardened concrete is 
related to R^/c and Ra/c as well as the curing conditions. For example, a high 
value of Ryj/c yields a low strength concrete owing to an increase in porosity 
that occurs because of the excess water present during hydration; consequently, 
most concrete is mixed with an initial water-to-cement ratio ranging from 0.30 
to 0.60. 

The bottom end of the concrete sample is immersed in water, where we 
impose the following Dirichlet boundary condition 

0(O,t)=0„,,-^^^, (15) 
Pg 

which states simply that the sample is fully saturated at x = 0. We also assume 
perfect sink conditions on the aqueous species, so that 

Cj(0,t) = forj = a,/3,g. (16) 

In typical experiments, the concrete sample is coated on the sides and top 
face with a sealant (such as epoxy) that prevents any transport into or out 
of the sample. This supports our ID approximation and allows us to impose 
homogeneous Neumann boundary conditions 

_(£,t)^0 and ^(i,t) = 0, (17) 

where j = a, (3, q. These conditions are equivalent to imposing a zero flux 
because the boundary condition on 6 at x = L requires that the convective flux 
component is zero. We note in closing that no boundary conditions are needed 
for Cg because it is governed by an ODE. 

3.4 Reaction rates 

The reaction terms are specified using notation introduced by Papadakis et al. 
wherein the rate of generation Vj [g/cm^ day] of species j = a,P is 




. = k,C, , (18) 
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with kj [day~^] a rate constant, Uj a power-law exponent, and C° the initial 
concentration (all given in Table [T]). The total rate of generation of C-S-H due 
to the alite and belite reactions [g / cm? day] is 

mesh ( ra_ rp \ 



where TOq, and mesh are molar masses of alite, belite and C-S-H respectively. 
A power-law reaction mechanism similar to (|18|) has also been employed in other 
models of cement chemistry [37l |42l [17] . 



3.5 Moisture diffusion coefficient 

Following the approach used for biofilms in ,15] we take the effective diffusivity 
to be a separable function of the form 

Di9,e)^ip^'^/^D*{d), (20) 

where the influence of porosity on clogging appears as a power law in the quan- 
tity If — ^o~_^g"" ■ Clement et al. initially assume that the ratio of porosities 
obeys ip = {R/'Ro)™, where R and Ro represent the corresponding pore radii 
and TO is an empirical constant. They then take two very common functional 
forms of the constitutive laws for porous media (namely, the van Genuchten and 
Brooks-Corey relationships) and show that the hydraulic conductivity in both 
cases satisfies K/Ko = (^(5»"+4)/2m. ^.j^g diffusivity must obey a similar relation- 
ship since it is proportional to K. By comparing with experimental data from a 
wide range of soils, Clement et al. find their power-law fit to be insensitive to the 
specific choice of to. They conclude that to = 3 is a reasonable approximation, 
which corresponds to the exponent 19/6 used above. 

The question remains whether these relationships applied successfully to 
biofilm growth in soils are also applicable to C-S-H gel formation in concrete. 
It is certainly true that the physics governing the two processes are very dif- 
ferent. Nonetheless, models for formation of C-S-H are based on the premise 
that the gel precipitates as outgrowths from the surface of cement grains [M] 
which is analogous to the manner in which biofilms accumulate on soil particles. 
Furthermore, the derivation above uses only spatially averaged quantities and 
hence makes no assumption about any specific pattern of biofilm growth. We 
therefore conclude that the 19/6 rule should also be applicable to concrete. 

Turning now to D* (9) , we observe that many other studies of water transport 
in concrete and related porous media [151 [321 [Mj approximate the diffusivity 
by an exponential function of saturation 

D*{0)=Ae^'^, (21) 

where parameters A [cm^ / day] and B are empirical constants. Lockington 
et al. |34| performed extensive experiments which showed that a number of 
building materials may be characterized by a universal exponent represented by 
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the rescaled parameter B = B [9 max — Omin) whose value lies between 4 and 6; 
other work [39] suggests that B could be as low as 2 and as high as 8. Note that 
these parameters lead to very rapid variations in diffusivity over the physical 
range of saturations (by at least three orders of magnitude) which distinguishes 
water transport in concrete from that of many other common porous media. 

3.6 Choice of base case parameters 

The numerical simulations in this paper focus on reproducing experimental re- 
sults reported by Barrita et al. [H [3] and specifically the concrete sample they 
refer to as "mixture 3." We begin by selecting a representative set of parameters 
for a "base case" simulation, but since not all of the required data is provided 
in these references we have had to estimate certain values using other literature 
sources. The parameters are summarized in Tableland we comment below on 
a number of the more critical choices: 

Sample geometry. We have taken the model domain to have length L = 
10 cm which is consistent with the cylindrical samples of concrete used in |3j. 

Water transport coefficients. The maximum water content is 9max = 0.067, 
which is equal to the initial porosity for the base case and is consistent with 
measured values reported in the concrete literature [Sj [31] . The residual water 
content is taken to be a small positive number because concrete is typically not 
totally free of water unless it has been artificially dried ,19] and in practice some 
small amount of water is typically trapped within the porous concrete matrix; 
specifically, we choose a value of 9„iin — 0.04 by estimating the minimum water 
content from plots in [4]. We take the diffusion parameter B — 100, which 
is chosen so that the rescaled quantity B = B{9max — 9min) = 2.66 (for the 
base case and other simulations performed later) is consistent with the range of 
values mentioned in Section [331 The value oi A — 0.0028 then follows by fitting 
the simulated wetting curves to Barrita's experimental results (more details are 
provided in Section 14. ip . 

Diffusion coefficients for aqueous species. Since the alite and belite actually 
dissociate and diffuse as ions, the best we can do is to use an approximation 
that represents the diffusivities in some averaged sense. We begin with the dif- 
fusivities of the ionic constituents Ca^"*", 0H~ and H2Si04~, which are equal to 
0.68, 4.6, and 0.43 cm? /d respectively [13], and compute an appropriately-scaled 
harmonic average of approximately 1.0 cm? jd for both alite and belite (follow- 
ing the development in [TB]). The diffusion of ions in cementitious materials is 
known to be reduced by a factor ranging from 10~^ to 10~^ [13] that depends 
on the pore structure and cement composition; in the absence of any better 
information we choose a factor of 10~^ after which = Djj = 0.01 cm? jd. 
The C-S-H gel does not diffuse in ionic form, and since no data is available in 
the literature regarding its diffusion coefficient we have chosen to simply take 
the same value Dq = 0.01 err? jd. This is not so much of a concern, since we 
investigate later on in Section [46l the effect of varying Dj and demonstrate that 
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Table 1: Parameter values corresponding to the base case. 



Symbol 


Description 


Value 


Units 


Reference 


Pw 


Liquid water density 


1.0 


•J ! 




Pn 


C-S-H eel density 


2.6 


q/cm^ 
•J 1 


Allen et al. [2] 


P CG7TL 


Cement density 


2.83 


gj cnt" 


Barrita et al. [3] 


Pagg 


Aggregate particle density 


2.6 


gjcrv? 




ma 


Alite molar mass 


228.3 


g/mol 




171(3 


Belite molar mass 


172.2 


g/mol 






Water molar mass 


18.0 


g/mol 




rricsh 


C-S-H molar mass 


342.4 


g/mol 
cm^ / day 
cm^ / day 
cm^ / day 
cm^ / day 




Da 


Alite difFusivity 


0.01 




Dp 


Belite difFusivity 


0.01 




D, 


C-S-H (aq) diffusivity 


0.01 




A 


Water diffusion coefficient 


0.0028 




B 


Water diffusion exponent 


100 


— 


Lockington et al. [31] 




Residual water content 


0.04 


— 


Barrita |4j 


Of 


Reaction cut-off 


0.04 


— 


Equal to d-^tn 




Alite reaction rate 


22.2 


day ^ 


Papadakis et al. [38] 


k(3 


Belite reaction rate 


3.04 


day^ 


ri 




Alite reaction exponent 


2.65 


— 


ri 


n(3 


Belite reaction exponent 


3.10 


— 


ri 


k'pVG.C 


C-S-H precipitation rate 


32.2 


day^^ 


Bentz [H] 




C-S-H dissolution rate 





day~^ 


ri 


V 


Water stoichiometry 


5 




Eqs. (dJ) and @ 


Rw/c 


Water-to-cement ratio 


0.333 




Barrita et al. [3] 


Ra/c 


Aggregate-to-cement ratio 


2.86 


- 


n 


LJa 


Alite mass fraction 


0.65 




ri 


Up 


Belite mass fraction 


0.17 




ri 




C3A mass fraction 


0.11 




ri 


u 


Alite hydration fraction 


0.60 




Tennis & Jennings [U] 


fp 


Belite hydration fraction 


0.20 




ri 


f. 


C3A hydration fraction 


0.72 




ri 


AVa 


Alite volume change 


0.233 


cm^ / g 


ri 




Belite volume change 


0.228 


cm? / g 


ri 


AVj 


C3A volume change 


0.555 


cm? / g 


ri 


L 


Length of sample 


10.0 


cm, 


Barrita et al. [3] 




Initial C-S-H (aq) concentration 





g/ cw? 






Initial C-S-H (gel) concentra- 
tion 





g/crv? 





Derived 



'parameters: 

Initial alite concentration 0.145 

Initial belite concentration 0.076 

Initial porosity 0.067 

Maximum water content 0.067 



g/cm^ Eq. ^ 
g/crr? " 

Eq. ([H 
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the solution is relatively insensitive to the values of diffusivity. 

C-S-H composition. C-S-H takes on a whole range of possible forms rep- 
resented by the general formula Cy^2i^z and so can only be considered in an 
averaged sense. We take rUcsh = 342.4 g/mol as a representative molar mass 
corresponding to y = 2: = 3, which is consistent with many other studies. There 
is a correspondingly wide range of gel densities reported in the literature, from 
1.85 gjcnt' at the lower end [55] up to 3.42 gjcm? [H]; we have chosen an 
intermediate value of pg — 2.6 g/cm^ which is justified by recent work on C-S-H 
microstructure . 

Cement composition. According to information provided in [3| on concrete 
mixture 3, the mass fractions of silicate constituents in the cement are oja = 0.65, 
ujp = 0.17 and uj-y ~ 0.11, while the aggregate- and water-to-cement ratios are 
Rajc — 2.86 and R^jc — 0.333. The cement mixture also contains 30% by 
weight of fly ash, which is a lower-density pozzolanic additive that improves the 
strength and workability of the resulting concrete. Based on densities of 3.15 
and 2.08 g/cm? for Portland cement and fly ash respectively, this translates 
into an overall cement density of pcem — 2.86 gj err? . All concrete samples were 
moist cured for 7 days which allows us to estimate /„ = 0.60, //j = 0.20 and 

= 0.72 from the plot of hydration fractions versus curing time given in [44] . 
Finally, the aggregates used in all mixtures are a combination of both fine and 
coarse quartz materials, and so we take pagg = 2.6 gjcnt' which is representative 
of the dry particle density of sand. 

Alitc and belite reaction rates. There is considerable variation in rate 
parameters reported in the literature owing partly to the fact that many ex- 
periments are performed not on cement samples but rather under idealized 
equilibrium conditions in which reactants are in solution. We have therefore 
chosen our parameters based on the data provided in [38j . who proposed the 
mechanism (jl8|) along with reaction exponents = 2.65 and np — 3.10; how- 
ever if we use their values of ka — 1.01 and /c/j = 0.138, then our model exhibits 
negligible clogging. But in fact, the reaction rate coefficients reported in the 
literature vary by several orders of magnitude [46l [121 HI] and so this ambiguity 
has led us to use the reaction rates as fitting parameters. Specifically, we take 
ka — 22.2 and kp — 3.04, which lie within the range of published values while 
also maintaining the same ratio of k^/kp used in [38 (more details on the fitting 
procedure are provided in Section 14. ip . 

Precipitation and dissolution rates. Bentz |12j developed a model that 
assumes a linear hydration rate law with rate constant ranging from 0.264 to 
1.464 day~^ depending on R^/c- We choose the precipitation rate as the upper 
end of their range, but again scale using the same factor as the other reaction 
rates to obtain kprec — 32.2. We also take kdiss — following Bentz and others 
who neglect C-S-H dissolution. 
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4 Numerical simulations 



The governing equations are discretized in space using a centered finite volume 
approach wherein the domain is divided into N uniform ceUs having width 
h = L/N and centered at Xi — {i — l/2)h for i = 1,2, . . . ,iV. The discrete 
solution components, for example Ci(t) « Ca{xi,t), are approximations of the 
average value of the solution within each cell. Using this notation, the discrete 
approximation of the alite equation ([5]) is 



d{OiCi) _ Da ( Ci+i ~ Ci Ci — Ci-i 

^^1+1/2^+1/2 ^ '"i-l/2C''i-l/2 ,„ 



(22) 



where the quantities Ci±i/2 are approximations of the solution at the left (— ) 
and right (+) cell edges for which we use an arithmetic average Ci±i/2 = 
(Ci + Ci±i)/2. The discrete velocity at cell edges is written using the centered 
difference approximation of Darcy's law 

"i-1/2 — -i^y^i-l/2i£,i-l/2) 7 • 



The same approach is used to discretize the remaining conservation equations 
®, (Uni) and pT|) . In all cases, the equations corresponding to boundary 
cells 1 = 1 and N involve "fictitious" solution values located at points xq — —h/2 
and XN+i = L -\- h/2 which lie one-half grid cell outside the domain. The 
boundary conditions are discretized using second-order differences or averages, 
and are used to eliminate these fictitious values in terms of interior solution 
components. 

The resulting semi-discretization is fully second order accurate in space and 
leads to a system of 5A^ ordinary differential equations for the discrete solution 
values which we then integrate in time using Matlab's stiff solver odel5s. For 
all simulations, we use N — 100 cells and set both relative and absolute error 
tolerances for odelBs to 10~^. The equations are integrated to time i = 28 days, 
which requires approximately 40 s of clock time on a Macintosh PowerBook with 
a 1.67 GHz PowerPC G4 processor. 



4.1 Base case with and without reactions 

We focus on developing comparisons with the experiments of Barrita [H [3] who 
studied wetting of concrete cylinders with both water and isopropanol. The 
latter solute is particularly useful in such a study because the silicate compounds 
in concrete do not react with isopropanol as they do with water, and so the 
isopropanol results may be used to calibrate the diffusion parameter A with 
experimental data. 

In the absence of reactions {k^ — kp — kprec — 0) there is no change in 
constituent concentrations and so the problem reduces to a single nonlinear dif- 
fusion equation for the water content. It is well known that for an exponential 



15 



difFusivity of the form (PT|) with large B and small A, the diffusion equation 
has a solution which forms a steep front that progresses into the sample with 
speed nearly proportional to the square root of time; consequently, a plot of the 
isopropanol wetting front location versus t^^"^ should be a straight line, as indi- 
cated by the experimental data of Barrita reproduced in Fig. [2] (square points). 
The experimental results for water uptake also exhibit a linear trend over the 
first 8-10 hours, which represents a period over which the reactions have not 
yet begun to take hold and no significant clogging has occurred. However, there 
is a noticeable difference between the initial slopes of the isopropanol and wa- 
ter data, which is most likely attributable to variations in the concrete samples 
used, or differences in the capillary pressure or other transport properties for 
the two liquids. Therefore, we have fit our model to the first 8-10 hours from 
the water experiment instead of using the isopropanol data. 

We proceed by setting B = 100 and varying A until the slope of the wetting 
front curve best approximates that of the experimental data. This fitting yields 
an estimate of A = 0.0028 cm^/ day which is consistent with values reported by 
other authors such as [T] . A plot of the computed wetting front location (without 
reactions) is displayed in Fig. [2] alongside the corresponding experimental data 
for comparison purposes. In this and all successive computations, the front 
location s{t) has been estimated by identifying the point x where the water 
content comes to within some tolerance of 9mm] that is, s(t) = min{a; : 9{x, t) < 

0rmn + 0.0005}. 



Comp., no reactions 

Comp., with reactions 

□ Expt., isopropanol 

o Expt., water 




Time [day 



Figure 2: Wetting front location s(t) for the base case computations both with 
reactions (solid line) and without (dashed line), which should be compared to 
the experimental data with water (circular points, taken from f^). The cor- 
responding experimental data for isopropanol (square points) are also included 
for comparison purposes. 

Moving now to the case of water uptake including hydration reactions, it 
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remains to choose an appropriate scaling of rate constants in order to best 
match the location of the stalled wetting front in experiments. As mentioned 
earlier in Section 13.61 we scale the reaction and absorption rates ka , kp and 
kprec simultaneously with the same value, while holding their ratio constant at 
1.01 : 0.138 : 1.464. This procedure yields the rates ka = 22.2, fc^ = 3.04 
and kprec = 32.2 for which the computed water content profile is displayed for 
comparison purposes in Fig. [3l We observe that incorporating the effects of 
hydration reactions and clogging due to C-S-H gel formation clearly causes the 
wetting front to stall a short distance inside the sample. 




Figure 3: Plots of computed water content for the base case parameters, both 
with and without reactions. The various solution profiles correspond to 10 
equally-spaced times over the 28 days of the simulation. 

Plots of concentrations and gel-modified porosity are provided in Fig. 21 
which indicate how transport of reactants into the sample is initially dominated 
by diffusion (for which the front propagates with velocity proportional to i^/^) 
but then later stalls as C-S-H forms and is precipitated near the lower end of the 
sample. The onset of clogging can be clearly seen in the gel concentration plots 
where Cg exhibits a peak slightly behind the stall location, while the porosity 
drops to its minimum value (approx. 0.13) within an interval containing the 
wetting front and Cg peak. It is worthwhile noting that diffusion and reaction 
processes continue to occur even after the front stalls - most noticeably ahead 
of the wetting front - owing to the presence of residual pore water, although 
this process continues at a much slower rate. We emphasize that although the 
capillary percolation threshold = 0min corresponds to the point where water 
can no longer move by capillary action, there is still sufficient water available 
for the aqueous components to diffuse (since 6 represents the physically bound 
or absorbed water as well as capillary water). 

The trends shown here suggest that onset of clogging occurs in the interior 
of the sample to the right of the inflow boundary. This effect can be attributed 
to a large initial influx of water at x = that dissolves the alite and belite 
near the boundary transporting them some distance downstream before the gel 
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precipitates. This result is consistent with [J] who reported high values of water 
flux within the first few hours of their experiments. 

We emphasize here that similar stalling behavior has been reported by sev- 
eral other authors performing experiments on porous building materials [43[ 1321 
[55] although these authors attributed this behaviour to an anomalous diffusion 
mechanism. Our primary aim here has been to show that a similar phenomenon 
can arise from pore clogging caused by hydration of residual silicates in concrete. 

The formation of a wetting front and subsequent stalling due to pore clog- 
ging are strongly dependent on two components of our model: the porosity 
dependence in the diffusion coefficient which drops to zero as e ^ dmin and the 
"shut-off" factor {9 — 6r)~^ appearing in the reaction terms. To illustrate the 
impact of omitting either effect, we present two additional simulations. First, 
if the porosity correction factor is removed from D{6,e) in Eq. ((20|) . then the 
wetting front propagates as if there were no clogging at all. This is evident by 
comparing the plot of water content in Fig. [5^ with that from the non-reactive 
case in Fig.|3^. There is clearly no visible effect on the front motion, even though 
a significant level of C-S-H gel builds up due to the reactions (see Fig. [SJd). 

To investigate the effect of slightly relaxing the cut-off factor {9 — 9r)^ in the 
reaction terms, we replace the zero cut-off with a small positive value of 5 x 10~^ 
when 9 < 9r- The wetting front still stalls as indicated in Fig. [5^; however, 
reactions occur over the entire domain giving rise to a significant concentration 
of C-S-H gel to the right of the front and a corresponding small reduction 
in saturation below This effect may be attributed to self-desiccation; 

however, with no more guidance in how to determine the value of the cut-off 
parameter, we leave the study of this effect as a possible avenue for future work 
and retain the factor {9 — 6r)^ as is. Taking a larger value of the cut-off (close 
to 9min in magnitude) can lead to runaway reactions and instabilities in the 
numerical method. 

We conclude from these last two simulations that in order for our model to 
give a reasonable picture of clogging observed in re-wetting experiments, there 
must be some retarding of liquid transport through a porosity dependence in 
the water diffusivity, and furthermore the reactions must include a shut-off term 
similar to {6 — 6r)^ , although a small positive reaction rate might be allowed 
near the residual saturation. 

4.2 Grid refinement study 

To ensure that our numerical simulations are computing a consistent solution 
that converges with the expected order of accuracy, we performed a grid refine- 
ment study. The base case simulation was repeated on successively finer grids 
with TV = 25, 50, 100, 200, 400, 800, and the solution on the finest grid is treated 
as the exact solution. The solution error was estimated using the discrete £2 
norm of the difference in aqueous C-S-H concentrations \\C^ — C^"'^''* 11^2 ! 
solution component would suffice, but we choose Cq because it often displays 
the greatest variations. The results are summarized in Table [21 and the ratio 
between successive errors indicates that the solution appears to be converging 
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0.03 


a 






0.01 ' ' ' ' ' ' ' 1 

1 2 3 4 5 6 7 
Height [cm] 

e. Gel modified porosity, e. 

Figure 4: The remaining base case solution profiles corresponding to Figs. [2] 
and[3)3. In each plot, the solution is displayed at 10 equally-spaced time intervals 
over 28 days. The arrows on each plot indicate the progression of curves in the 
direction of increasing time. 
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3 4 5 
Height [cm] 

a. Water content. 




b. C-S-H gel concentration. 



Figure 5: Solution with no porosity dependence in the diffusion coefBcient, 
exhibiting an absence of clogging. 



0.07 

0.065 

0.06 

i 0.055 

c 

o 

S 0.05 
0.045 
0.04 




1 2 3 4 5 
Height [cm] 

a. Water content. 




1 2 3 4 5 6 7 
Height cm 

b. C-S-H gel concentration. 

Figure 6: Solution computed by replacing {9 — 9r)^ « ma.x{9 — 6*^,5 x 10~^), 
corresponding to the situation when reactions do not entirely shut off at the 
residual saturation. 
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at a rate that is at least second order, as expected. 

Table 2: Grid refinement study. The order is calculated as log2 (ratio). 



No. of points (N) 


£2-eTTOT 


Ratio 


Order 


25 


0.019 


2.12 


1.08 


50 


0.0087 


4.30 


2.10 


100 


0.0020 


5.56 


2.48 


200 


0.00036 


6.27 


2.65 


400 


0.00058 







4.3 Sensitivity to alite/belite reaction rates 

In this section we vary the reaction rate parameters and kp to investigate 
the effect of changes in the individual rates as well as the relative importance 
of the two reaction routes leading to production of C-S-H gel. To this end we 
hold k/^ constant and scale ka by the factors 0, 0.1 and 10, and then repeat 
the same procedure for k/s. The resulting solutions are displayed in Figs. [7] 
and [H from which we see that the clogging seen in the final solution is very 
sensitive to changes in both rates. The results in both cases are similar, with 
effect of alite being more pronounced; this is not surprising considering that the 
initial concentration of alite is significantly larger than that of belite (refer to 
values of C° and in Table [T|). We also note that if the alite reaction rate is 
taken small enough, then no stalling occurs and the wetting front propagates 
essentially unhindered into the sample; the same is not true of the belite rate 
since there is still enough alite being hydrated to cause significant clogging. 
The sensitivity to reaction rates demonstrated by these results points to the 
importance of obtaining accurate estimates of the rate parameters. 

4.4 Sensitivity to precipitation rate 

Since there is some uncertainty in the choice of the precipitation rate, it is help- 
ful to consider the effect of changes in kprec- We ran three additional simulations 
with kprec = 0.0, 3.22 and 322 and compared those to the base case in Figure [HI 
The kprec = case is identical to the case displayed in Fig [3^ (without reac- 
tions) and from the remaining results it is clear that the solution is relatively 
sensitive to the choice of precipitation rate. We have done our best to choose 
a value of kprec consistent with C-S-H precipitation rates in the literature, but 
there is potentially much to be learned by taking a more detailed look at the 
precipitation process and including more details about this and other reaction 
mechanisms in the model equations. 

4.5 Sensitivity to dissolution rate 

We have so far assumed that the formation of C-S-H (gel) is an irreversible pro- 
cess and no dissolution occurs, which is consistent with assumptions made in 
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1 2 3 4 5 6 7 1 2 3 4 5 

Height [cm] Time [day ] 

a. Final water content. b. Wetting front position. 



Figure 7: Water content and wetting front location for different values of the 
alite reaction rate, fc^. In this and all succeeding figures, the base case is plotted 
using a solid black line and highlighted in the legend using "*" . 




Heigiit [cm] Time [day ] 

a. Final water content. b. Wetting front position. 



Figure 8: Water content and wetting front location for different values of the 
belite reaction rate, kp. 
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Figure 9: Water content and wetting front location for different values of the 
C-S-H precipitation rate, fcprec- 



many other models. Since our focus is on the phenomenon of re- wetting wherein 
time scales are much longer than typically considered for initial hydration reac- 
tions, it is helpful to consider the effect of incorporating a non-zero dissolution 
rate constant kdiss- To this end, we considered values of kdiss — 1 and and 10 
day~^ and compared the resulting solutions in Fig. [TOl which clearly indicates 
that only for the largest value of kdiss is there any appreciable effect on the 
wetting front position, although the water content does show some deviations 
at smaller values of kdiss ■ These results support our assumption that dissolution 
has a negligible effect on the solution when kdiss *C kprec- 




Height [cm] Time [day""^] 

a. Final water content. b. Wetting front position. 



Figure 10: Water content and wetting front location for different values of the 
dissolution rate, kdiss- 
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4.6 Sensitivity to constituent difFusivity 

We next investigate the effect of changing the diffusion coefficients for the aque- 
ous alite, belite and C-S-H species. We note that our model ignores transport 
and reaction of individual ionic species and instead approximates the diffusive 
transport by employing an effective diffusion coefficient for each constituent 
which may not be entirely representative of how the individual ions would move 
in response to concentration gradients in solution. Fig. [TT] demonstrates that 
changes in the diffusion coefficient by several orders of magnitude have some 
effect on the steepness of the wetting front and the distribution of constituents 
behind it, but have very little influence on the location of the front itself. 




a. Final water content. b. Wetting front position. 

Figure 11: Water content and wetting front location obtained by varying the 
diffusivities Da, Dp and Dq. In each case depicted, all three diffusivities are 
scaled by the same constant factor. 



4.7 Sensitivity to aggregate density 

The aggregate materials typically used in concrete include sand and gravel of 
varying coarseness, all of which have different density. In practice, a combination 
of various aggregates is frequently used and so we next investigate the effect of 
variations in the aggregate density. Fig. [12] compares the solution when pagg is 
varied between 2.4 and 2.8, and shows that even such seemingly small changes 
in aggregate density can have a measurable effect on clogging; in particular, 
as Pagg increases, the degree of clogging experienced decreases. We therefore 
conclude that an inaccurate value of the aggregate density parameter could 
lead to incorrect results. 

4.8 Effect of changes in cement mixture 

Most concrete is mixed with a water-to-cement ratio Rw/c lying somewhere 
between 0.3 and 0.6. It is well known that when Rm/c is too large the resulting 
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a. Final water content. b. Wetting front position. 



Figure 12: Water content and wetting front location for different values of the 
aggregate density, pagg. 



concrete can be weak and so a smaller Ry^jc is desirable in general. On the other 
hand, if there is too little water then the cement can become unworkable or there 
may even be insufficient pore water to fully hydrate the silicates in the hydration 
process. Consequently, optimizing concrete strength and durability requires a 
fine tuning of the initial water content. We have simulated the effect of changes 
in composition by taking parameters as listed in Table [3l which correspond to 
mixtures numbered 1 through 4 from [3]. Outside of the variations in R^/c and 
Ra/ci a major difference between the various mixtures is the presence of fly ash 
(in mixtures 2 and 3) or silica fume (in mixture 4). Both of these low-density 
cement additives have the effect of reducing the value of Pcem, which in the case 
of mixtures 2 and 3 can change the resulting porosity e° significantly. 

Table 3: Composition of cement mixtures taken from Tab. 2], with computed 
results compared in Fig. [T21 



Mixture 


Pceni 


■^w / C 


Ra/c 


e° 


1 


3.15 


0.599 


5.39 


0.113 


2 


2.62 


0.364 


3.13 


0.074 


3 (base) 


2.83 


0.333 


2.86 


0.066 


4 


3.07 


0.297 


3.12 


0.045 



The resulting numerical solutions are compared in Fig. [T3] from which it is 
clear that the initial porosity (as determined by the concrete mixture) can have 
a major impact on water transport. We note in particular that mixture 1 (with 
the largest value of e°) exhibits no clogging, while the low value of e° in mixture 
4 leads to very limited water transport, with the wetting front stalling much 
closer to a; = 0. Indeed, Barrita et al. [3] observed in experiments that their 
mixture 4 exhibited a much earlier onset of clogging than the other concrete 
samples, an effect that is clearly captured in our simulations. However, there 
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remains some discrepancy in that experiments on mixture 1 exhibited a stalled 
wetting front, while our simulations show no clogging in this case. 




Heighl [cm] 

a. Final water content. 




b. Wetting front position. 



Figure 13: Water content and wetting front location obtained for various cement 
compositions, using mixtures 1-4 in [3]. 



5 Conclusions and future work 

We have developed a model for the transport and reaction of water and other 
reactant species in hardened concrete subject to re-wetting. Numerical simula- 
tions support our hypothesis that hydration of residual silicates and subsequent 
formation of C-S-H gel may be responsible for the clogging phenomenon ob- 
served in experiments, which is the main contribution of this paper. 

We investigated the sensitivity of the solution to changes in a number of 
model parameters, from which we can conclude that the reaction rate parame- 
ters (specifically k^, fc^ and kprec) have the most impact on the solution. These 
are precisely the parameters which are most difficult to ascertain owing to dis- 
crepancies in the published literature, and in particular the lack of values for 
reaction rates in actual concrete as opposed to idealized values obtained for 
silicates prepared in aqueous solutions. Consequently, more work is required to 
ensure that inputs to our model are consistent with actual concrete re-wetting 
scenarios. 

In addition to obtaining better estimates of the model parameters, there are 
a number of extensions to the current model which may significantly improve its 
predictive power. We expect that the greatest impact may be had by replacing 
the simple precipitation process embodied in our rate parameter kprec with a 
more realistic reaction mechanism that takes into account details of the C-S-H 
microstructure and hydration which have recently been uncovered. Possible 
examples include: 

• Incorporating the dynamics of individual ionic species through the addi- 
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tion of new transport equations and reaction kinetics along the lines of [5T] 
or [32]. 

• Investigating the hypothesis put forward in [461 that hydration kinetics 
is a two-stage process, consisting of an early accelerated hydration step 
followed by a slower hydration reaction that dominates in the longer term. 
They suggest that this two-stage kinetics might arise from effects of either 
C-S-H microstructure or precipitation kinetics, cither of which could be 
considered in detail by appropriate modifications of our model. 

• Separating the C-S-H gel into two forms characterized by different densi- 
ties as suggested in 43', '44' , where the lower-density gel is thought to be 
primarily responsible for changes in porous structure. Taylor et al. [43j 
also mention the importance of swelling in the cement matrix during initial 
cement hydration, which is an effect we have so far neglected. 

• Chemical shrinkage and the associated phenomenon of self-desiccation, 
which are known to have a significant impact on initial cement hydra- 
tion [40] . 

It may prove useful to incorporate other aspects of porous transport that are 
commonly seen in modelling studies of ground water aquifers or oil reservoirs, 
but have yet to be applied to the study of concrete. For example, capillary 
hysteresis has been identified as an important aspect of cement hydration [S] 
and results from the soil sciences community [371 130] could certainly be applied 
in this context. The issues raised in [211 surrounding the impact of variable 
porosity on models of multi-phase transport should also be applicable to cement 
and concrete. Our model can be easily adapted to study other stages in the life 
of concrete such as initial hydration, carbonation, aging or degradation. Finally, 
the approach we have developed here would also be applicable to the study of 
other transport phenomena such as polymer flooding in enhanced oil recovery, 
where chemical reactions and solution-dependent parameters are important. 

Acknowledgements 

This work was supported by grants from the Natural Sciences and Engineering 
Research Council of Canada and the MITACS Network of Centres of Excellence. 
JMS was supported by a Research Fellowship from the Alexander von Humboldt 
Foundation during a visit to the Fraunhofer Institut Techno- und Wirtschafts- 
mathematik in Kaiserslautern. We thank Dr. Jesus Cano Barrita (Instituto 
Politecnico Nacional, Oaxaca, Mexico) for many insightful discussions and for 
providing the experimental data. We are also sincerely grateful to the four 
anonymous referees whose extensive comments have significantly improved this 
work. 



27 



References 



[1] H. Akita, T. Fujiwara, and Y. Ozaka. A practical procedure for the analysis 
of moisture transfer within concrete due to drying. Mag. Concr. Res., 
49(179):129 137, 1997. 

[2] Andrew J. Allen, Jeffrey J. Thomas, and Hamlin M. Jennings. Composition 
and density of nanoscale calcium-silicate- hydrate in cement. Nature Mater., 
6:311-316, April 2007. 

[3] P. F. de J. Cano Barrita, B. J. Balcom, T. W. Brcmncr, M. B. MacMillan, 
and W. S. Langley. Moisture distribution in drying ordinary and high 
performance concrete cured in a simulated hot dry climate. Mater. Struct., 
37:522-531, 2004. 

[4] Prisciliano Felipe dc Jcsiis Cano Barrita. Curing of high-performance con- 
crete in hot dry climates studied using magnetic resonance imaging. PhD 
thesis, University of New Brunswick, Fredericton, NB, November 2002. 

[5] B. Bary and A. Sellier. Coupled moisture-carbon dioxide-calcium transfer 
model for carbonation of concrete. Cement Cone. Res., 34:1859-1872, 2004. 

[6] Z. P. Bazant and L. J. Najjar. Drying of concrete as a nonlinear diffusion 
problem. Cement Cone. Res., 1:461-473, 1971. 

[7] Jacob Bear. Dynamics of Fluids in Porous Media. Dover, New York, 1988. 

[8] James J. Beaudoin. Why engineers need materials science. Concr. Inter- 
nat, 21(8):86-89, August 1999. 

[9] Arnon Bcntur. Ccmcntitious materials - Nine millennia and a new century: 
Past, present, and future. J. Mater. Civ. Eng., 14(l):2-22, 2002. 

[10] D. P. Bentz. Effects of cement PSD on porosity percolation and self- 
desiccation. In B. Persson and G. Fagerlund, editors. Self- Desiccation and 
Its Importance in Concrete Technology II, pages 127-134. Lund University, 
1999. 

[11] Dale P. Bentz. A three-dimensional cement hydration and microstructure 
program. I. Hydration rate, heat of hydration and chemical shrinkage. Re- 
port No. NISTIR 5756, Building and Fire Research Laboratory, National 
Institute of Standards and Technology, Gaithersburg, MD, 1995. 

[12] Dale P. Bentz. Influence of watcr-to-ccmcnt ratio on hydration kinet- 
ics: Simple models based on spatial considerations. Cement Cone. Res., 
36(2):238-244, 2006. 

[13] J. Billingham, D. T. I. Francis, A. C. King, and A. M. Harrisson. A multi- 
phase model for the early stages of the hydration of retarded oilwell cement. 
J. Eng. Math., 53:99-112, 2005. 



28 



[14] J. D. Birchall, A. J. Howard, and J. E. Bailey. On the hydration of Portland 
cement. Proc. Roy. Soc. Land. A, 360:445-453, 1978. 

[15] Michael Chapwanya, Stephen B. G. O'Brien, and J. F. Williams. A ID 
bio-clogging model in a phreatic aquifer. In preparation, 2008. 

[16] T. P. Clement, B. S. Hooker, and R. S. Skecn. Macroscopic models for 
predicting changes in satm-atcd porous media properties cause by microbial 
growth. Ground Water, 34(5): 934-942, 1996. 

[17] M. M. Y. Delmi, A. Ait-Mokhtar, and O. Amiri. ModeUing the coupled 
evolution of hydration and porosity of cement-based materials. Constr. 
Build. Mater., 20:504-514, 2006. 

[18] Daniels Farrington and Robert A. Alberty. Physical Chemistry. Wiley, 
1966. 

[19] Daniele Ferretti and Zdenek P. Bazant. Stability of ancient masonry tow- 
ers: Moisture diffusion, carbonation and size effect. Cement Cone. Res., 
36:1379-1388, 2006. 

[20] E. J. Garboczi and D. P. Bentz. Computer simulation of the diffusivity of 
cement-based materials. J. Mater. Sci, 27:2083-2092, 1992. 

[21] W. G. Gray and C. T. Miller. Examination of Darcy's law for flow in porous 
media with variable porosity. Env. Sci. Tech., 38:5895-5901, 2004. 

[22] M. W. Grutzeck. A new model for the formation of calcium silicate hydrate 
(C-S-H). Mater. Res. Innov., 3:160-170, 1999. 

[23] C. Hall, W. D. Hoff, S. C. Taylor, M. A. Wilson, Beom-Gi Yoon, H.-W. 
Reinhardt, M. Sosoro, P. Meredith, and A. M. Donald. Water anomaly in 
capillary liquid absorption by cement-based materials. J. Mater. Sci. Lett., 
14:1178-1181, 1995. 

[24] Christopher Hall. Anomalous diffusion in unsaturated flow: Fact or fiction? 

Cement Cone. Res., 37:378-385, 2007. 

[25] Christopher Hall, W. D. Hoff, and M. Skeldon. The sorptivity of brick: 
dependence on the initial water content. J. Phys. D, 16:1875-1880, 1983. 

[26] Torben C. Hansen. Physical structure of hardened cement paste. A classical 
approach. Mater. Constr., 19(114):423-436, 1986. 

[27] Daniel Hillel. Environmental Soil Physics. Academic Press, 1998. 

[28] Hamlin M. Jennings, Jeffrey J. Thomas, Julia S. Gevrenov, Georgios Con- 
stantinides, and Franz- Josef Ulm. A multi-technique investigation of the 
nanoporosity of cement paste. Cement Cone. Res., 37:329-336, 2007. 



29 



[29] Jcspcr Kildsgaard and Peter Engesgaard. Numerical analysis of biological 
clogging in two-dimensional sand box experiments. J. Contam. HydroL, 
50:261-285, 2001. 

[30] J. B. Kool and J. C. Parker. Development and evaluation of closed- form 
expressions for hysteretic soil hydraulic properties. Water Resour. Res., 

23(1):104 114, 1987. 

[31] Rakcsh Kumar and B. Bhattacharjcc. Porosity, pore size distribution and 
in situ strength of concrete. Cem,ent Cone. Res., 33:155 164, 2003. 

[32] Michel Kiintz and Paul Lavallee. Experimental evidence and theoretical 
analysis of anomalous diffusion during water infiltration in porous building 
materials. J. Phys. D, 34:2547-2554, 2001. 

[33] Frederick M. Lea. The Chemistry of Cement and Concrete. Edwin Arnold, 
Glasgow, third edition, 1970. 

[34] D. Lockington, J.-Y. Parlange, and P. Dux. Sorptivity and the estimation 
of water penetration into unsaturated concrete. Mater. Struct., 32:342-347, 
1999. 

[35] D. A. Lockington and J.-Y. Parlange. Anomalous water absorption in 
porous materials. J. Phys. D, 36:760-767, 2003. 

[36] Marc Mainguy, Claire Tognazzi, Jean-Michel Torrcnti, and Frederic Ade- 
not. Modelling of leaching in pure cement paste and mortar. Cement Cone. 
Res., 30:83-90, 2000. 

[37] S. A. Meier, M. A. Peter, A. Muntean, and M. Bohm. Dynamics of the 
internal reaction layer arising during carbonation of concrete. Chem. Eng. 
Sci., 62:1125-1137, 2007. 

[38] V. G. Papadakis, C. G. Vayenas, and M. N. Fardis. A reaction engineering 
approach to the problem of concrete carbonation. AIChE J., 35(10):1639- 
1650, 1989. 

[39] Leendert Pel. Moisture transport in porous building materials. PhD thesis, 
Technische Universiteit Eindhoven, February 7, 1995. 

[40] Bertil Persson. Self-desiccation and its importance in concrete technology. 
Mater. Struct, 30:293-305, 1997. 

[41] S. J. Preece, J. Billingham, and A. C. King. On the initial stages of cement 
hydration. J. Eng. Math., 40:43- 58, 2001. 

[42] Anna V. Saetta, Bernhard A. Schrefler, and Renato V. Vitaliani. 2-D model 
for carbonation and moisture/heat flow in porous materials. Cement Cone. 
Res., 25(8):1703-1712, 1995. 



30 



[43] S. C. Taylor, W. D. HoflF, M. A. Wilson, and K. M. Green. Anomalous water 
transport properties of Portland and blended cement-based materials. J. 
Mater. Sci. Lett, 18(23):1925-1927, 1999. 

[44] Paul D. Tennis and Hamlin M. Jennings. A model for two types of calcium 
silicate hydrate in the microstructure of Portland cement pastes. Cement 
Cone. Res., 30:855-863, 2000. 

[45] Jeffrey J. Thomas and Hamlin M. Jennings. Effects of d20 and mixing on 
the early hydration kinetics of tricalcium silicate. Chem. Mater., 11:1907- 
1914, 1999. 

[46] F. Tzschichholz and H. Zanni. Global hydration kinetics of tricalcium 
silicate cement. Phys. Rev. E, 64:016115, 2001. 

[47] WBCSD. The cement sustainability initiative: Our agenda for ac- 
tion. World Business Council for Sustainable Development, Available at 
|http : //www ■ wbcsdcement . org , July 2002. 



31 




Height [cm] 



